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Mott transitions are studied in the two-dimensional Hubbard model by a non-perturbative theory of correlator 
projection that systematically includes spatial correlations into the dynamical mean-field theory. A nonzero 
second-neighbor transfer yields a first-order Mott transition at finite temperatures ending at a critical point. 
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The Mott transition has been one of the fun- 
damental issues in strongly correlated electrons 
pp. In spite of substantial progress |1I2| . its full 
understanding in low dimensions remains open. 

Here, we propose a non-perturbative theoret- 
ical framework of correlator projection method 
(CPM) 131. The CPM combines the operator pro- 
jection method ^ with the dynamical mean-field 
approximation (DMFA) by restoring spatial cor- 
relations along the coarse-graining concept in the 
energy-momentum space with increasing the or- 
der of projection. It does not suffer from a cluster- 
size problem. Then, we obtain the Mott transi- 
tion phase diagram of the two-dimensional (2D) 
Hubbard model in the parameter space of the 
Coulomb repulsion U ^ the second-neighbor trans- 
fer t' , the temperature T (scaled by the nearest- 
neighbor transfer t) and the filling (n). A first- 
order Mott tansition surface with a critical end 
curve appears at half filling. Doping into the 
Mott insulator yields a diverging compressibility 
and/or a tendency to phase separation. 

We consider the 2D Hubbard model. The for- 
malism of CPM [3j is as follows: The operator 
projection leads to series of Dyson equations, 

E„_i(.^,k) = G„_i(c.,k)4"^"-^^ (1) 
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G„-i(^,k) = l/[^-4"'"^-S„(c.,k)] (2) 

with n ~ 1, 2, • • •. Go = So is the electron 
Green's function and Ei is the self-energy part 
with £^^'°^ = 1. Other coefficients up to the sec- 
ond order n — 2 are given by the Hartree-Fock 
dispersion e'^'^^ = — tk — /i — C/(n}/2, a source 
of Hubbard band splitting e^'^^ = U'^M with 

M = (n)(2 - (n})/4, and 4^'^^ = - A* - 
[/(I — {n)/2) — Ecov/M. Here, /i is the chemical 
potential determined from the filling {n). Scor = 
- I]x',s*x,x'((l/2 - nx -s)(4,sCx',s + h.c.))/2 is 
a local energy shift due to a correlated hopping 
process with rix.s = s^x.s where Cx.s and c]^ ,. 
are the electron creation and annihilation oper- 
ators at a site x with a spin s, respectively. 
and are the Fourier transforms of the trans- 
fer ix,x' and tx,x' = ^x,x'(('^x'«x')/4-H (SW ■ S^') - 
(AjjAx'))/-^, respectively, where rix, 5^ and Ax 
represent the charge, spin and local-pair opera- 
tors, respectively, it introduces k dependences of 
I]i(ci;,k) mainly through the superexchange pro- 
cess and produces antiferromagnetic precursors. 
Here, we determine it from Ref. 0. 

To reproduce a Mott insulator in particle- 
hole asymmetric cases, the local dynamics of 
S]2(w,k) must be obtained in a sufficiently cor- 
rect manner. Then, instead of decoupling ap- 
proximations |4I6| ■ we adopt the following gen- 
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eralized DMFA: (i) We calculate a local normal- 
ized self-energy part GijocC^) = j^J^k^ii^^^) 
from an arbitrary E2,ioc(w) by replacing E2(w,k) 
with E2joc(w) in Eq. (2), (ii) calculate the 
Weiss self-energy partSi{ijj) = e'^'^'^'' /[G^\^^{lo)-\- 
S2,ioc('^)], and (iii) obtain the Weiss Green's 
function Qa{uj^]<i) = — e^^'^"* — Siiuj)]. (iv) 
Within the iterative perturbation scheme, we cal- 
culate a new I]2joc(w) as a Fourier transform of 

'^2,\oc(j ) = lv3 I]k,k',q rk,k',q^/o('r , ^+^G{){t , k'- 

q)C'o(-T,k') where Tk.k'.q = 4[tk' -^k'-q-ik+q^- 
^k + 4eeor/((ri)(2 - {n)))f / {{n){2 - (n))). This 
loop continues to (i) until convergence is reached. 

This second-order CPM reproduces a Mott in- 
sulator at low T's for i' = at half filling and its 
single-particle spectra and momentum distribu- 
tions show a remarkable similarity to numerical 
results 17]. As we introduce t' ^ with U and 
{n) ~ 1 being fixed, the system undergoes a first- 
order Mott metal-insulator transition (MIT) at a 
low T accompanied by a discontinuity in the dou- 
ble occupancy (n^ui). In the space of {U,t',T), 
this first-order MIT boundary forms a surface . 
Within the surface with U being fixed, the jump 
in (n^-ni) decreases with increasing T and disap- 
pears above a critical temperature Tcr; the first- 
order surface ends at a critical curve at T > 0. 
Then, without any spontaneous symmetry break- 
ing as in two dimensions at T > 0, the Mott 
transition shows a similarity to a liquid-gas phase 
transition, as in Refs. |8l2j . The MIT boundary 
at T agrees with numerical results 

Doping carriers into Mott insulators yields 
a metallic state characterized by an enhanced 
charge compressibility k = d{n) /dfj,. The present 
results of /i — (5 curves near half filling are shown in 
Fig. ^with 6 = l — {n). The slope dfi/d{n) repre- 
sents 1/k. Then, from high T's, k increases with 
decreasing T and shows a divergence towards the 
filling-controll MIT at low T's. Even an instabil- 
ity to a phase separation appears as emergence 
of K < at lower T's. However, the issue if the 
phase separation is an artifact of the present ap- 
proximation requires a further intensive study. 

In summary, using the CPM that systemati- 
cally includes spatial correlations into the DMFA, 
the MIT phase diagram of the 2D Hubbard model 



has been clarified in the ([/, space with a 

first-order MIT surface with a T > critical end 
curve for t' ^ at half filling as well as a diverg- 
ing charge compressibility and/or a tendency to 
a phase separation. 




Figure 1. fi-S curves obtained with U/t ~ 4. The 
bold circles lie in the Mott-insulating side of the 
MIT surface at half filling. 
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